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We study the Hamiltonian approach to Yang-Mills theory in Coulomb gauge in Rayleigh- 
Schrodinger perturbation theory. The static gluon and ghost propagator as well as the potential 
between static colour sources are calculated to one-loop order. Furthermore, the one-loop /3- function 
is calculated from both the ghost-gluon vertex and the static potential and found to agree with the 
result of covariant perturbation theory. 

Q\ ', PACS numbers: ll.lO.Ef, 12.38.Bx 

■ Keywords: Coulomb gauge, Hamiltonian approach, perturbation theory 

o ■ 
(N : 

■ I. INTRODUCTION 

^: 

' In recent years, there has been a renewed interest in Yang-Mills theory in Coulomb gauge both in the continuum 
psj ! 0, H, 0, 0, H, H S Q and on the lattice [1, [13, [HI ■ This gauge has several advantages over the frequently used Landau 
Ps| ' gauge. Among these are: the use of the physical degrees of freedom (at least in QED) and the explicit emergence of 
] a static colour charge potential. Furthermore, in this gauge the form factor of the ghost propagator represents the 
' ^ ' ■ dielectric function of the Yang-Mills vacuum The disadvantage of this gauge is, of course, that it is non-covariant, 
' which is considered a drawback for perturbation theory, and renormalisability is yet to be proven. Recently, there has 
Q^r. been much activity in the Hamiltonian approach to Yang-Mills theory in Coulomb gauge [l], 0, H, [l^ [3, [l^ IJ^ . 
By means of a physically motivated ansatz for the vacuum wave functional, a variational solution of the Yang-Mills 
Schrodinger equation has been accomplished 0, 0, 0] ■ The so-called gap equation resulting from the minimization 
of the vacuum energy density was converted into a set of Dyson-Schwinger equations for the static gluon, ghost and 
I [ Coulomb propagators. Due to the particular ansatz for the vacuum wave functional used, which grasps the essential 
i infrared physics, the resulting gluon propagator does not yield the proper ultraviolet (UV) asymptotics known from 
perturbation theory. The reason is that the three-gluon vertex does not contribute to the vacuum energy density for 
On 1 the ansatz of the wave functional considered. 

\l ■ In the present paper, we carry out a thorough perturbative calculation in the Hamiltonian approach to Yang-Mills 
, theory. We will use the standard Rayleigh-Schrodinger perturbation theory and carry out the calculations up to 

■ and including one-loop order. In particular, we will produce the correct one-loop expression for the perturbative 
' /3-function which has been found in the functional integral approach ^, i8j] before. Because of asymptotic freedom, 

we thus have a correct description of the UV regime (to one- loop precision). With the correct perturbative results at 
hand, we show how the variational approach has to be modified to yield the correct UV asymptotics. Furthermore, the 
perturbative calculations performed in the present paper can be extended to calculate systematic corrections to the 
non-perturbative variational vacuum solution and to construct a complete basis for the Hilbert space of Yang-Mills 
^ , theory, which is required for the calculation of the partition function or free energy. In a forthcoming paper, the 
' latter will be minimized to achieve a description of the deconfinement phase transition at finite temperatures in the 
Hamiltonian approach. 
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II. PERTURBATIVE EXPANSION OF THE YANG-MILLS HAMILTONIAN 



A. The Yang— Mills Hamiltonian in Coulomb gauge 

In the absence of matter fields, the Yang-Mills Hamiltonian in Coulomb gauge reads [l3| 



YM 



(x) + i F^^ (x) F^. (x) 



(1) 



where x is a vector in d space dimensions, is the transverse gauge field operator, 11° is the transverse momentum 
operator satisfying the canonical commutation relations 



and 



(2) 



(3) 



is the field strength tensor {g is the coupling constant, /"'"^ are the structure constants of the su(A^c) algebra and 
j^ic _ jabcjfb gauge field operator in the adjoint colour representation). Furthermore, 



J[A] = Det(-£> • d)/Dei{~d'^) , bf = 5^^ d, + gAf 
is the Faddeev-Popov determinant of Coulomb gauge, and 

i^'^^[A](x,y) = U-b ■ i-d^) {-b ■ d)-' 



ah 
x,y 



(4) 



(5) 



the so-called Coulomb kernel. 

The Coulomb gauge diA^ = is implemented in the scalar product of the Hilbert space of Yang-Mills wave functionals 
by the Faddeev-Popov method, yielding for the matrix elements of observables 0[A,n] 

(6) 



(7) 



(^1 1 1 V^2> - jvA J[A] ri [A] 0[A, n] V2 [A] , 
where the integration is over the transverse gauge fields only. 

It is convenient to remove the Jacobian from the integration measure by defining 



The transformed Hamiltonian H is then obtained from Eq. ([T]) by replacing the momentum operator H by the 
transformed one 



n^(x) = = n^(x) + ^^1^ 



(8) 



This yields 



YM 



iHf(x)n,"(x) + ^F^.(x)^^5(x) 



^ 1 d'^xd'^yir(x)nf(x)F'^''[A](x,y)if(y)H^^(y). 



(9) 



In the following, Dirac's bra-ket notation will refer to the transformed space where no Faddeev-Popov determinant 
occurs in the functional integration measure. 

For perturbation theory, it will be convenient to expand the gauge field in Fourier modes. We use the following 
conventions 



(10) 



where the transformed fields satisfy the canonical commutations relations 



(11) 



and where tij (k) 



kikj is the transverse projector in momentimr space. To simplify the notation, we 



introduce the following shortcuts 



dp ': 



{2nY 



(12) 



B. Expansion of the Hamiltonian 

Expanding the Hamiltonian Hyu © hi powers of the coupling constant g, thereby using — 1 + 0{g^) and 



we obtain 



where the unperturbed Hamiltonian 



H = Ho+gH,+g^H2+Oig''), 



(H^(x)) -^^(x)aM^(x) 



(13) 
(14) 

(15) 



is the Hamiltonian of QED except for the extra colour index of the gauge field. The first order term g Hi arises from 
the expansion of the magnetic energy J d'^xFf^ and is given by the three-gluon vertex 

Hl^J^ jaia.a3 j ^fc^ ^(^1 + ks + kg) T(l, 2, 3) ^1(1) A{2) A(3) . (16) 

Here we have introduced the shorthand notation A{\) = A°^(ki) and T(l,2,3) carries the (totally antisymmetric) 
Lorentz structure of the three-gluon vertex, 



T(l,2,3) := tiij(ki)ti2/(k2)ti3m(k3)[(5j7(fc2 - ^l)m + ^im(fc3 - ^2)^ + ^im(fcl " ^3);] 

Finally, the second-order term 



(17) 



H2 = - y dfci . . . dfc4 <f (ki H- k2 + k3 + k4) 

^-^l^^p^ A{1) A{2) A{-i) A(A) ' ^'^'-^'^'^ 



(ki-f k2)^ 



^(1)H(2) A(3)H(4) 



(18) 



C. 



contains besides the usual four-gluon vertex (first term in the bracket) also a contribution from the Coulomb term 

(second term in the bracket) arising from the expansion of the Coulomb kernel ([5|) . Note that the Coulomb term is 

already 0{g'^), see Eq. Q, so that to the order considered we can replace the Coulomb kernel ([5]) simply by its bare 

form (— 9^)~^. The last term in Eq. (fT5|) is an irrelevant constant arising from the expansion of the second term in 

Eq. (fT5|) . Since such a constant does not influence the wave functional, we will skip it in the following. 

Since we will use dimensional regularisation, in order to preserve the dimension of the dressing functions we will 

replace 



(3-d)/2 



(19) 



with pi being an arbitrary mass scale. 



C. The unperturbed basis 



The perturbative vacuum state V'o[^] = (^|0) is given by 

{A\0) =AAcxp|-i J dfcA^(k)t,, (k)|k|A^"(-k)| 



(20) 



Up to the colour index of the gauge field, this is precisely the exact vacuum wave functional for QED without fcrmions. 
This state is the lowest energy eigenstate of Hq (fTS|) with energy 



Eo = {Nl - 1 



d-l 



&k Ikl 



(21) 



To restrict the Coulomb gauge field to its transverse degrees of freedom, it is convenient to introduce the eigenvectors 
of the transverse projector ty (k) in d spatial dimensions with eigenvalue 1 



t.,(k)cj(k)=cf(k), 



(22) 



where a labels the d—1 different eigenvectors. From kitij{\<i) = it follows immediately that these eigenvectors are 
orthogonal to k 



Assuming the normalisation 



fc,<(k) = 0. 



cr(k)c[(k)=J. 



these vectors satisfy the "completeness" relation in the transverse subspace 

cnk)cf (k)=t,,(k). 



(23) 
(24) 
(25) 



Since the transverse projector tij (k) is a real symmetric matrix, it can be diagonalised by an orthogonal transformation 
and hence the eigenvectors c°'(k) can be chosen real, which we will assume below. Let us mention, however, that in 
d = 3 the cf (k) are usually chosen as the circular polarisation vectors, which are complex. 
The transverse components of the gauge field and their momenta are then given by 



A^(k) := <(k)A^(k) , n^(k) := <(k)n^(k) . 
In view of Eq. (jlip they satisfy the commutation relation 

[A^(k),n^(p)] =z5.x^<f(k + p). 

In the transverse components of the gauge field, the perturbative vacuum state ([201) reads 

(A|0) =AAexp|-i j dfcA^(k)|k|A^(-k) 
This state is annihilated by the operator 

a^(k) = yS[A^(k) + ^n^(k) 
which together with its hermitean conjugate 

«r(k) = 



<(-k)-— n^(-k) 



fulfills the usual Bose commutation relation 

[a^(k),a^t(p)] ^^abj^^5(k_p) 



, a^(k)|0)=0, 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



The unperturbed Hamiltonian Eq. (|15p is diagonalised by the transformations (P^ . ([50]) 



Ho = Eo+J dk\k\a''J{k)aU\^), 
and accordingly the eigcnfunctions of i/o ^^'^ the muhiple giuon states defined by 



\9i92 ...gN)^ allH^i) Cj(k2) . . . aZH^N) |0) 



with energies 



AT 



En 



(32) 



(33) 



(34) 



i=l 



There are two ways to proceed now: Given the unperturbed basis p3l) . one can express the gauge field A and its 
momentum operator 11 in the perturbations, Eqs. ()16p and (|18p . in terms of the creation and annihilation operators 
aj., flo-. Alternatively, since the perturbation Hi,2 is expressed in terms of the gauge field and its momentum, one may 
wish to express the unperturbed basis states through A and 11 using Eq. ((30|) . This yields 



l5> 



^^^(-k)|0). 



for a one-particle state, and 
151,52) 



,5(^1 



|0), 



(35) 



(36) 



for a two-particle state, where the second term on the right hand side of the last equation arises from the canonical 
commutation relations. States with more gluons have similar additional contraction terms. These contact terms 
ensure that the free n-gluon states are orthogonal to each other and to the unperturbed vacuum. It turns out that 
the contact terms simply eliminate from the matrix elements of observables the contractions of gauge field operators 
stemming exclusively from the wave functionals. We can therefore use the simplified representation 



I5i, 



,9n) = 



|0), 



(37) 



with the additional calculational rule that in the evaluation of matrix elements of the form (0| 0[A, 11] |gi, . . . , g„), the 
A fields in the states Eq. ([57|l must not be contracted with each other but only with A, H occurring in the observable 
0[A,U] . 



D. Expansion of the vacuum wave functional 



In Rayleigh-Schrodinger perturbation theory, the leading order (in the coupling constant) corrections to the vacuum 
wave functional are given by 



10)^ 



^ 1 (gi....g„|gi|0) 



E 



n! |ki| + ... 

1 1 
^ |ki| + ... 



|k„| 



(51... 5«|i?2|0) + (<?!... 5n|i^l|0) |5l ■ • ■ 5n 



(38a) 
(38b) 



where Hi and H2 are defined in Eqs. (|16p. (fT5| . and the factors 1/n! avoid multiple counting due to identical gluons 
(summation over colour and polarisation indices and integration over the momenta is implicit). 

It is now straightforward to calculate the perturbative corrections to the vacuum wave functional. Since Hi contains 
three field operators and is antisymmetric in both colour and Lorentz indices, in view of Eq. ([57| it is clear that only 
three-gluon states will contribute to |0)' Eq. (|38ap . yielding 



5(ki -f k2 -f ks) 



r(l,2,3)v4(-l)A(-2)A(-3)|0). 



(39) 



As a consequence, the second term in |0)' ' , Eq. (|38b|) . receives contributions from up to six-gluon states. Furthermore, 

— (2) 

since H2 contains terms with up to four field operators, the first term in |0) will receive contributions from two- and 
four-gluon states. However, it turns out that up to order only two-gluon states contribute to the static propagators, 
so that 



Nr. — I dfci dfc2 ^ , 2) + Fc{l, 2) + 2)] A{-1) A{-2) |0) 



10) ... 3 

where we have introduced the abbreviations 



Fb(1, 2)^-JSq 



^c(l,2)^/d.^ 
Fi(l,2) = -i J dk 



(kl-q) = 



T(-l,3,4)r(2,3,4)5(k2+k3 + k4) 



|ki| + |k3| + |k4| 



Iksl |k4| 



(40) 

(41a) 
(41b) 
(41c) 



Here, Fb{1, 2) arises from the four-gluon vertex, Fc{l, 2) from the Coulomb term, and Fi(l, 2) from the second-order 
contribution of the three-gluon vertex, see the last term in Eq. (|38bp . We have also used the shorthand notation 

In the Raylcigh-Schrodinger perturbation theory used in Eqs. psp . the correction to a wave function of a given order 
is chosen to be orthogonal to the unperturbed wave function. 



(0|0) =0,i>l 



(42) 



As a consequence, the wave functions obtained in a given order of perturbation theory are not properly normalised. 
We are interested in the vacuum wave functional up to second order 



\n)^N-n \0)+g\or' +g'\0f' +Oig') 



(43) 



Calculating the normalisation constant Afn up to order making use of Eq. ()42p , we find for the properly normalised 
vacuum wave functional up to this order 



\n) 



\0)+g\0)' +g'\oy +Oig') 



(44) 



In the following sections we will use this perturbative expansion of the vacuum wave functional to calculate various 
static propagators. 



III. GHOST PROPAGATOR 



The Green's function (or inverse) G[A] of the Faddeev-Popov operator is defined by 

- Df[A]i^)df G''^[^](x,y) = <5"^ <5(x - y) . 
Expanding G[A] in a power series in the coupling constant 

G'^^[A](x,y) = ^5"Gf[A](x,y) 



71=0 



we get from (j45|) the recursion relation 



dl Gf+i[^](x,y) = irWar GJl''[^](x,y) , n > 



together with the initial condition 



5,^GS''[A](x,y) = <5''''<5(x-y), 



(45) 
(46) 

(47a) 
(47b) 



which defines the unperturbed static ghost propagator G'g''[A](x, y) as the Green's function of the Laplacian, which 
in momentum space reads 



GS''[A](k) = 5"^Go(k) = ^ 



(48) 



It is diagonal in colour space and independent of the gauge field. With this property, it follows from Eq. (|47ap that 
each term G„[^] contains a product of n gauge field operators. 

The ghost propagator is defined as the expectation value of the inverse Faddeev-Popov operator in the vacuum state 



G{x,y) = {n\G[A]{^,y)\n) . 



(49) 



Since G[A] does not depend on the momentum operator 11, it does not change under the transformation to the "radial" 
Hilbert space, i.e. G[A] = G[A], see Eq. ([T]). Furthermore, contrary to the Faddeev-Popov operator, the vacuum 
expectation value of its inverse, the static ghost propagator, is translationally invariant. 

Inserting the expansions (|46|) for G[A] and (I44|) for jfi) into Eq. (I49|) it is possible to show that many terms vanish or 
cancel, so that the ghost propagator reduces to 



G(x, y) = Go (x, y) + g^0\ G^ [A] (x, y)\0) +0{g'). 
The second-order term in Eq. ([SO)) can be evaluated by means of the recursion relation (|47ap . yielding 



Gf[A](x,y)- Go (A- d) Go (A- d) Go , Go = (-9^) 



ab 



(50) 



(51) 



The vacuum expectation value of Eq. (j5ip can be expressed through the bare static gluon propagator, which in view 
of Eq. (HOI) reads 



|A(1) A(2)|0) = 5°i°^5(ki 



2|ki| 



Then for the ghost form factor Dc(k) defined by 

G'(x,y) 

we get the following expression at one-loop order 

D,{k) = l + .gV' 



£'c(k) 



2 ,,3-d ^<l_ ^ fc.; kj tjj (q) 



2k^ 



(k-q)2 |q| 



(52) 



(53) 



(54) 



The integral ([51)1 is standard and can be evaluated in dimensional regularisation with ci = 3 — 2e in the usual way, 
yielding 



Nr 



(47r) 



4 /I 



8 



28 



- - - In — - 7 - - ln2 + — + Oie 



(55) 



IV. GLUON PROPAGATOR 



A. Gluon propagator in the Hamiltonian approach 



The full static gluon propagator is defined in momentum space by 

s{k + p) Di'{k) = {n\An\^)A%p)\n). 

With the expansion (|44p for the vacuum functional, the nonvanishing terms up to order 0{g^) are 



'>|A(l)A(2)|0r + (01^(1) A(2)|0r + '10|A(1)A(2)|0) 



(56) 



(57) 



FIG. 1; Diagrammatic representation of Eq. (|59|l . The curly and the fuU lines represent, respectively, the static gluon and 
Coulomb propagators. 



where Dq is the tree- level propagator given in Eq. (l52|) . As for the static ghost propagator, there are no terms of 
0{g), since the connected pieces of {0\AA\0} ' are given by the expectation value of five gauge field operators, which 
vanishes. 

The normalisation factor [second term in the first line of Eq. ()57p ] cancels the disconnected piece of the first term in 
the second line. From Eq. (|57p it is also clear that the contributions to |0)* ' with more than two-gluon states do not 
contribute to D{k) to the order considered. 

The matrix elements in Eq. ((S7)) can be straightforwardly evaluated using the results of Sec. [TTl Moreover, the last 
term in the second line of Eq. (IF7|) can be expressed in terms of the preceeding one since 

>|A(1)A(2)|0) = >|At(-l)At(-2)|0) = ((0|A(-l)A(-2)|0f 

For the gluon form factor -DA(k) defined by 



2|k| 



(58) 



we get at one-loop level 



N, 



8(<i- I)k^ 



d<7 



- 3rf + 3 - (k • q)2 



Nr. 



4(d~I) 



ty(k) iy(q) 



(k-q)2 



1 



(59) 



8((i- l)k^ 



dgdp 



T(p,k,q)(27r)'^<5(p + k + q) 



(IpI 



|q|j|pMq| 



with k = k/|k| and 



T(ki,k2,k3) =tr|r(l,2,3)|2, 
where T(l, 2, 3) is defined in Eq. P7|) and the trace is taken in Lorentz space, yielding 

T(k, q, p) = 2(p2 + q2 + k^) [d - 2 - (p • q)(q • k)(p • k)] + 

+ 2(d - 1) [(p • q)(p • k) + (q • p)(q • k) -f (k • q)(k • p) - p • q - q • k - p • k] ^ 
- 2{d - 3) [(p • q)2(p2 + q2) + (q . k)2(q2 + k^) + (p . k)2(p2 + ^2^] 



(60) 



(61) 



The diagrammatic representation of Eq. (j59p is shown in Fig.[TJ^ 

The first integral in Eq. (l5^ is a tadpole term, which vanishes identically in dimensional regularisation. Moreover, in 
the last line of Eq. (I59p we can carry out one of the momentum integrals due to the (5-function, yielding 



Nr. 



l + /.(k) 



4((i- l)k^ 

Nr 

2(d- l)k^ 



d - 2 -f (k ■ q) 
(k-q)2 



S(k,q) 2|kH-|q| + |k-q| 



q||k-q| (|k| + |q| + |k-q|)^ 



(62) 



^ Notice that these graphs are not standard Feynman diagrams. 



with 



q) = Uiik) tj„i{q) tfc„(k - q) [Sijkk - SjkQi - Sikkj] [Simkn - S„inqi - Snik^] 

{d-2)k^ q2 + (k-q)2" 



l-(k.q: 



(rf-l)(k^ + q^) 



(k 



(63) 



The first integral in Eq. (|62p . i.e. the contribution /c(k) from the Coulomb kernel, can be evaluated in d = 3 
dimensions by means of standard techniques, yielding 



2s 



Ic{k) = 



(47r)2- 



4 
15 



-In- 



15 



In 2 



188 
225 



Oie) 



(64) 



Unfortunately, the gluon loop /g(k) [second integral of Eq. (j62p ] is highly non-trivial. It would probably be possible 
to evaluate /g(k) using partial differential equations techniques similar to the ones used in Refs. 0, Sj- Instead of 
using these techniques, we will show that the integrals in Eq. ([5^ are the same as the ones treated in Refs. 0, @|i ^^'^ 
we will use the result of those papers. 



B. Static gluon propagator from the Lagrangian approach 

In the Lagrangian-based functional integral approach to Yang-Mills theory in Coulomb gauge considered in Refs. 0)11], 
the full (energy-dependent) propagator has the form 

{Anp)A%k)) = {27rf+U{p + k)6-'u,{k)W{k4,k), (65) 

where VF(fc4,k) can be expressed (in Euclidean space) as 

W{k4, k) = k2j^]^2 ■ (66) 

Here we have introduced the dressing function £'yi^(fc4,k), which measures the deviation of the propagator from the 
tree-level form. We are interested here in the static or equal-time propagator, the quantity considered in the Hamil- 
tonian approach, which is obtained from W{ki, k) by integrating out the temporal component of the 4-momentum 

W{k)^ J^W{k,,k). (67) 

At tree-level (where Daa = 1) this yields 

Woik) ^ ^ , (68) 

which is precisely the static tree- level gluon propagator of the Hamiltonian approach, see Eq. ((52|) . For sake of 
comparison with the Hamiltonian approach, Eqs. (j58p and (j62p . we also express the full equal-time gluon propagator 
((H7)) by a dressing function 5^^(k) 

Wik)=:DM^)Woik)m£^ (69) 

The two dressing functions (form factors) in and ([H^ are related by 

Z)..(k).2|k|/f^%i^. (70) 
J fc| + k 

With the dressing function of the energy-dependent propagator DAAik4,k) given in Refs. 0, [1], we can calculate the 
dressing function DAAik) of the equal-time propagator. We will now show that the so obtained equal-time dressing 
function of the gluon propagator Eq. ([70)1 coincides, at one-loop level, with the gluon form factor of the Hamiltonian 
approach, defined in Eq. 



In Refs. 0, [1| the dressing function DAA{k4,i(.) is evaluated at one-loop level, with the result-^ 



DAAih-.'ii-) 



2N, 



d'^qdqi E(k,q) 



{d-1) J {2TTY+i fc2 g2 (fc _ q)2 



d'^qdq4 Uj(k)Uj{(i) 2n 
l_K4 — q j 



(71) 



{d-1) J {2TTY+i fc2 g2 (k - q)2 



Here S is the kernel obtained by contracting the three-gluon vertex, defined in Eq. (j^]) . Inserting Eq. ((7T|) in Eq. ([70]) 
and also performing the loop integration over q^ one obtains for the equal-time dressing function Daa(^) (j70p precisely 
the form factor DA^k) of the static gluon propagator of the Hamiltonian approach, Eq. (I62p . We have thus shown 
that the Lagrangian-based functional integral approach yields the same equal-time gluon propagator as the time- 
independent Hamiltonian approach, at least to the order considered. We have checked that this equivalence does also 
hold for the (HH) correlator, but it does not hold for the (AH) correlator. The reason is that in the time-dependent 
Lagrangian approach the (AH) correlator is odd under time reversal and thus the corresponding equal-time correlator 
vanishes, while in the Hamiltonian approach the static (AH) correlator is constrained by the canonical commutation 
relation not to vanish. 

In Refs. the loop corrections to the gluon dressing function Eq. ([55|) were calculated in dimensional regularisation. 

With d — 3 ~ 2e the result reads 



DAA(A4,k) = 1 



(47r)2-s 



1 

- - 7 - In ^ 



+ 4 m 



1 

- - 1 

z 



llz-3z" 



- ln(l - 



64 , 



1 

2^ 



14 

y 



3 



Ois) 



(72) 



where z = fc|/k^ and the functions f{z), g{z) are defined by 



/(z) -41n2 



arctan y/z 



Injl + zt) 
iz Jo Vt{l + zt) 

l{z) = 2 In 2 - ln(l + z) 



(73a) 
(73b) 



Using this result, after integrating Eq. 
static gluon form factor 



72|) over ^4 with the appropriate tree- level factor, see Eq. (ffOl) . we find for the 



DaO^) = 1 



(4^) 



2-e 



In- 



(74) 



where the ellipsis contains the finite constant terms. 

Finally, we remark that the ghost propagator in the Lagrangian-based functional integral approach 0, [1] defined there 
as the correlator (cc) for explicitly introduced ghost and antighost fields, is given by a function Wc{k4, k) in analogy 
with Eq. (|65p for the gluon propagator, which is independent of the temporal component k^, T/Fc(fc4,k) = H/c(k). It is 
easily seen, by integrating over the temporal component of the loop momentum, that M^c(k) coincides with our result 
dSll) for G(k) = ^(k)^^. 



V. THE GHOST-GLUON VERTEX AND THE /3-FUNCTION 

In the Hamiltonian approach the ghost-gluon vertex is given by 

{n\An^)G''^'[A]{y„y,)\n) =: 

J d^zi d'^Z2 d-^za A, (x, zi) G(yi, Z2) G(y2, Z3) rf =(zi; Z2, Z3) , ^'^^^ 

where D is the full gluon propagator of Sec. IIV Al and G the static ghost propagator as defined in Sec. IIIII Eq. ([75]) 
tells us to calculate the vacuum expectation value (AG[A]) and then "cut off" the external legs, in order to extract 



The results of and 0, being evaluated in the, respectively, first and second order formalism, are at first sight not identical. However, 
it is not difficult to show that they indeed agree. 




FIG. 2: Diagrammatic representation of the one-loop correction to tiie ghost-gluon vertex. The left and right picture correspond 
to, respectively, the first and second integral in Eq. (I77|l . The arrow shows the flow of the loop momentum I. 



the irreducible component. In Eq. (j75p we already used the fact that ghost and giuon propagators are colour diagonal 
to every order in perturbation theory. 

The lowest-order contribution to Eq. (j75|) is the bare vertex r*-"-*, given in momentum space by 



rf)'^''^(k;p,q) = -*5r''^p., 



(76) 



(all momenta defined as incoming). With arguments similar to the ones we used in the evaluation of the ghost and 
gluon propagators, it is not difficult to show that the terms contributing to the next non-vanishing order are 



\AG2{A\\^) 



|AG2[A]|0) -^(0|AG3[A]|0)iPi, 



where the subscript IPI means that only the irreducible terms have to be considered. The evaluation of these matrix 
elements is straightforward and yields 



rp-(k;p,q).-.5V^^^f {/d^^ 



{£-p) 

i Pj Qm tjm (^) 

£)2(q + £)2 " 
r,j„i(k, p-.£,q + 



i^\p~£\\ci + e\ |k| + |p-.g| + |q + .6| 



(77) 



where T is the Lorentz structure of the three-gluon vertex (fT7|) . The diagrammatic representation of Eq. ([77|l is given 
in Fig. [21 The integrals in Eq. ([77]) are UV finite in c? = 3 spatial dimensions, hence they are independent of the 
scale when evaluated at a symmetry point (after factorising the momentum that carries the Lorentz index i). Then 
5^ DaDI is a renormalisation group invariant, at least to the present order. From Eqs. ([55]) and ([Ti]) we obtain in 
Coulomb gauge 

n 2 



1 



9 



(47r)2-e 

Nr. 



1 

- + 

e 

II 1 



I 



N, 4 1 
(47r)2-'^ 3 e 



9 



(47r)2-^ 3 e 

Consequently, must have a 1/e pole with coefficient (— lI/3)A^c, and since, at the one-loop level, the coefficient of 
the 1/e pole of is the (gauge invariant) first coefficient /3o of the /3-function, 

dg 



/3(5) = 



1 



(9 In /I (47r) 



■/3o 5-^ + 0(5^), 



we find 



which is the correct value. 



/3o = -y A^c, 



VI. THE POTENTIAL FOR STATIC SOURCES 



(78) 



(79) 



Until now, we have considered pure Yang-Mills theory without external colour charges. If static external sources are 
included, the Hamilton operator becomes 



H — Hym 
. 9' 



2 



I d'xd'y p:;,(x)F'^''[A](x,y)Af(y)n?(y)-HAr(x)n,^T(x)F'^''[A](x,y)p:;,(y) 



(80) 



where Hym is the Hamihonian ^ of the Yang-Mills sector considered in the previous sections and p°„(x) is the 
external charge density. Since 

^ ^ ' " ' ^ ' ^ 2 SA- (x) ' * ^ ' ~ ' ^ > 2 5 A- (x) ' 

the 5\nJ^ /5A contributions cancel in the terms linear in the external charge density and we can omit the tildes there, 
obtaining 



H^Hym + ^J d^xd^y p'^^{^)F'^''[A]i^,y)pl{y) + 

+ d''xd'*2/{p:^,(x)i^'^^[A](x,y),if(y)n,^(y)}. 



(81) 



({,} denotes the anticommutator.) We are interested here in a perturbative calculation of the static potential, i.e. 
the potential between static charges. For this purpose, we will follow the approach of Refs. [19, 20, 21], treating 
the external sources as perturbations to the pure Yang-Mills sector. Suppose we know the exact spectrum of Hym, 
i.e. its eigenstates \^n) and the corresponding eigenvalues -E^^. Second-order perturbation theory in the external 
sources then yields for the potential 

V-i..y) ^ (*.|F-|Al(x,y)|*„, - £ ^ (^.1 A^M /.'(y) |.,) ,,,, 

where we have introduced the quantity 

K'^i^) j d^'z {i^''^[A](x,z),i^(z)n^(z)}. (83) 

Since we are interested in the potential to order O(g^), we can replace all quantities in the second term in Eq. ([82|l 
by their unperturbed expressions. This yields 

^ ($o|j^°(x)|<l'Ar)(<&Ar|/v''(y)|$o) | (0| j^pin |gi , .92) ^ , 4, 

"T^^ E^-^, |k,| + |k,| (^4) 



where the unperturbed part of F'^''[A](x, y), Eq. in d = 3 spatial dimensions 

gab 

47r|x — y 



Fr(x,y)=rMMr^].,,-73^, (85) 



is the familiar Coulomb interaction. 

In the first term of Eq. ((5^ (the so-called colour Coulomb potential) we use the perturbative expansion of the vacuum 
wave functional Eq. and also expand the Coulomb kernel F'^^[A]{k, y), Eq. in powers of 5. Thereby many 
terms can be shown to vanish or cancel using similar arguments as in the evaluation of the static ghost propagator. 
In the end one obtains for the Coulomb potential 

{n\F^''[A]i^,y)\n) ^ FS''{^,y) + {0\Ff[A]{^,y)\0) +Oig') , (86) 

where F2''''[A](x, y) = 3Gf [A](x,y) is the 0{g^) part of the Coulomb kernel F[A], Eq. ©. 
For the potential dressing function v(k) defined by 

F'^''(x,y) = y'dfce*''(-'')^^^, (87) 

we obtain the following expression 

v[\.)-l + 9P ak^y ^^k-q)2|q| + 

_„2„3~d^ /•.^ t,,(q)t,,-(k-q)|k-q|-|q| (88) 
2k^J'"' |q| |k-q| + |q| 

= 1 + 'yi(k) + wii(k) . 



(a) (b) 



FIG. 3: Diagrams contributing to the static potential in second order perturbation theory, (a) and (b) represent, repectively, 
the second term of Eq. (f5^ and Eq. . 



The diagrammatic representation of this equation is given in Fig. [31 The first integral of Eq. (I88p comes from the 
expectation value of F2 [second term in Eq. (|86p]. while the second integral arises from the second-order term [Eq. ([HI])]. 
While the first integral in Eq. (|88p is standard, the second one is of the same type we encountered in the calculation 
of the static gluon propagator. Again, we will use the results of the Lagrangian-based functional integral approach. 
As argued by Zwanziger [22] , the potential between static colour charges should be related to the correlation function 
of the Aq field. For the corresponding dressing function defined by 



{Al{p)Al{k))^{27:t+U{p + k)& 
the result at one-loop order reads 0] 



ab Da{h,k) 



(89) 



d'^q (k'-2k-q)t,j(q)t,j(k-q) 
{27rr [k-q] [A:| + (lk-ql + lql)2] 



(90) 



where we have already performed the integration over the temporal component 54 of the loop momentum. Noticing 
that 



k^ - 2k • q = (k 

we see that at one-loop level 

Taking then the result for Dcr obtained in 7] 



- q)2 - q2 = (Ik - q] + 
v{k) = D,{ki^O,k) 



|q|) 



(91) 
(92) 



D^{ki,k) = 1 
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+ (3z ~ l)g{z) + -(1 + z)(l + 3z)/(z) + 0{e) 



(93) 



where z — k\/\s? and the functions f{z), g{z) are defined in Eqs. ([75]) . we get for the dressing function of the static 
potential 



i;(k) = 1 



Nr. 



(47r 



l2-e 



11 
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- - 7 - In — 
e /i^ 



(94) 



A few remarks are in order here. The equivalence between the static potential evaluated in the Hamiltonian approach 
and the propagator of the Aq field in the Lagrangian-based formalism is far from trivial, since in the Hamiltonian 
approach the Weyl gauge = is also imposed. Furthermore, note that the Hamiltonian static potential is not 
obtained as the equal-time component, but rather as the integral over the relative time of the {AqAq) correlator. 
Such an integral captures the instantaneous interaction through the first correction term in Eq. ([90]) as much as the 
retarded interaction via vacuum polarisation given by the second correction term there. 

Given the fact that the physical potential is g^V'^'^ , with F"'' given by Eq. ([57|) . t;(k) in Eq. ([55)) represents the form 
factor of the running coupling. From Eq. ([M]) we then find /3o — — llA^c/3, which is again the correct coefficient. 
Since this agrees with £'a(p)_D^(p), we find the same running coupling from the ghost-gluon vertex and the Coulomb 
potential, at least to the order considered. 



VII. RELATION WITH THE VARIATIONAL APPROACH 



In the variational approach considered in Refs. 0, £^ generahsation of Eq. (|20p was taken for the vacuum functional, 



{A\uj) ^AAexp^-- / ^A^(k)t,,(k)^(k)AK-k)^, (95) 



1 /• d^fc 

and the kernel w(k) was determined by minimisation of the vacuum energy density. In this non-perturbative case a 
complete basis can also be defined by Eqs. 150]) and ([331) with |k| replaced by ti-'(k). Considering the three-gluon vertex 
as a perturbation on top of the non-perturbative vacuum (|95|) . the energy functional gets an additional contribution 



3! t^(ki) + cj(k2) + tj(k3) 

Functional differentiation of this expression with respect to w(k) yields an additional term to the gap equation 

a;2(k)=k^+x'(k)+/°+/.(k)+ 

2 A^, [ d3g d3p r(k,q,p) (27r)3j(p + q + k) 2w(k) + tj(q) + t^(p) (97) 



^87 (2^)3 (2^)3 [c^(k)+a;(q)+u;(p)]c^(p)a;(q) w(k) + w(q) + Lj(p) ' 

where T(k, q, p) is given in Eq. (j^U)) and the integral terms x(k), and /cj(k) are defined in Ref. For k ^ cx3 
this additional term reduces to the integral /g(k), Eq. (|62p . found in perturbation theory, showing that Eq. ([W]) does 
indeed provide the right contribution to the variational form of the gap equation to produce the correct UV asymptotic 
behaviour, at the same time leaving the ghost dominated infrared sector untouched. The numerical solution of this 
modified gap equation is in progress. 



VIII. SUMMARY AND CONCLUSIONS 



In this paper we have studied the Hamiltonian approach to Yang-Mills theory in Coulomb gauge in Rayleigh- 
Schrodinger perturbation theory. The static gluon and ghost propagators as well as the potential between static 
colour sources have been calculated to one-loop order using dimensional regularization. The one-loop /3-function 
was calculated from the ghost-gluon vertex as well as from the static potential. In both cases the result known 
from covariant perturbation theory was reproduced. The unperturbed basis constructed from the eigenstates of 
the unperturbed Hamiltonian, which up to the colour index of the gauge field coincides with the Hamiltonian for 
QED (without fermions), was generalized to multiple quasi-gluon states on top of the non-perturbative vacuum wave 
functional used in the variational approach [3, IB]- Treating the three-gluon vertex, which is not captured by the 
variational wave functional used so far, as a perturbation on top of the non-perturbative vacuum, a modified gap 
equation was derived, which yields the correct (perturbative) ultraviolet asymptotics for the static propagators known 
from perturbation theory while, at the same time, leaving the (non-perturbative) infrared behaviour of the propagator 
unchanged. The multiple quasi-gluon basis constructed on top of the non-perturbative vacuum wave functional will 
also serve as a basis for calculating the partition function of Yang-Mills theory in the Hamiltonian approach in 
Coulomb gauge and investigating the deconfinement phase transition at finite temperatures. 
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